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1. Introduction 

Center projection in direct maximal center gauge is a procedure used to locate center 
vortices in lattice gauge field configurations. Results obtained from this procedure, 
which tend to support the center vortex theory of confinement, have been reported in 
many articles over the past several years. Last year, however, an alarming negative 
result was reported by Bornyakov, Komarov, and Polikarpov (BKP) in ref . [1] . Recall 
that direct maximal center gauge fixing is defined by the gauge transformation g 
which maximizes the sum 

R = ^Tr A [«U^x)] (1.1) 

where Tr^ indicates the trace in the adjoint representation, and 9 U fl (x) is the gauge- 
transformed lattice configuration. This prescription is simply Landau gauge fixing 
in the adjoint representation. There is no known method for finding the global 
maximum of R, but most previous work has used the over-relaxation scheme of ref. [2] 
to find local maxima (Gribov copies). The string tension that can be attributed to 
vortices, which are located using this scheme, agrees quite well with the asymptotic 
string tension of the full lattice configuration. This agreement is known as "center 
dominance," and is crucial to the argument that center vortices account for the 
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entire asymptotic string tension. In their work last year, BKP used an improved 
center gauge-fixing scheme based on simulated annealing, and obtained Gribov copies 
with consistently larger values of R than those obtained by over-relaxation method. 
The discouraging outcome of this improved gauge-fixing procedure is that center 
dominance is much less evident: the center projected SU(2) string tension at f3 — 2.5 
is some 30% lower than than the full asymptotic string tension. 

In this article we explain the origin of the difficulty found by BKP, and then 
propose and test a method which overcomes that difficulty. The realization that 
direct maximal center gauge must fail in the continuum limit is actually due to 
Engelhardt and Reinhardt in ref. [3], whose work predates the numerical results 
of BKP in ref. [1]. We have previously elaborated on the Engelhardt-Reinhardt 
reasoning in ref. [4]; the salient points will be reviewed in section 2 below. In brief, 
these authors pointed out that direct maximal center gauge can be understood as a 
best fit to a given lattice gauge field by a thin vortex configuration. Since the field 
strength of a thin vortex is highly singular, the fit must fail badly near the middle 
of the vortex. In the continuum limit, the bad fit to the gauge field in the vortex 
interior would overwhelm the good fit in the exterior region. As a result, in the 
continuum limit, the center projection obtained from a global maximum of R may 
reveal no vortices at all. 

Having diagnosed the illness, a cure is proposed in section 3. We work throughout 
with the SU(2) gauge group, although the analysis generalizes readily to SU(N). The 
new method calls for maximizing 

R M = J2 Tr[M T (x)U Af ,(x)M(x + /})] (1.2) 

x,n 

where Ua^{x) denotes the link variable in the adjoint representation, and where, for 
SU(2) gauge theory, M(x) is a real-valued 3x3 matrix which satisfies the unitarity 
constraint only as an average; i.e. 

^EE M 5(^ M ^) = ^ (i.3) 

V x j 

with V the lattice volume. The motivation for relaxing the local unitarity constraint 
is explained below. The real- valued matrix field M(x) which arrives at a global max- 
imum of Rm can be determined uniquely, in terms of the three lowest eigenfunctions 
of a lattice Laplacian operator. To determine the corresponding gauge transfor- 
mation, we relax M(x) to the closest SO(3) matrix-valued field Qa( x ) satisfying a 
related Laplacian condition. The SO (3) field is mapped to an SU(2) matrix- valued 
field g(x), which is used to gauge-transform the original lattice configuration. Center 
projection then determines the vortex locations. 

In the end, direct Laplacian center gauge is simply lattice Landau gauge, in the 
adjoint representation, with a particular choice of Gribov copy (different from the 
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global maximum of R). We will motivate this choice in sections 2 and 3. Numerical 
results obtained from the new procedure, regarding center dominance, vortex density 
scaling, precocious linearity, and the correlation of vortex locations with the values 
of Wilson loops, are reported in section 4. 

Direct Laplacian center gauge is closely related to the original Laplacian Landau 
gauge of Vink and Wiese [5], appropriately generalized to the adjoint representation. 
The main difference is in the mapping from M(x) to g(x). Direct Laplacian gauge also 
has much in common with the Laplacian center gauge, introduced by Alexandrou 
et al. [6] and de Forcrand and Pepe [7], in that both gauges require solving the 
lattice Laplacian eigenvalue problem in the adjoint representation. But the gauges 
themselves are rather different; they are motivated by different considerations (best 
fit in one case, gauge-fixing ambiguity in the other), and the numerical results are 
not the same. These points are further discussed in section 5. Section 6 contains our 
conclusions. 



2. Maximal Center Gauge and its Discontents 



We begin from the insight that lattice Landau gauge fixing is precisely equivalent 
to finding the best fit, to a given lattice gauge field U^(x), by a classical vacuum 
configuration g(x)g^(x + fi). The best fit vacuum configuration minimizes the mean 
square distance on the group manifold between vacuum and gauge field link variables 



d F = ^E Tr ( u i*( x ) ~ 9{x)g\x + fa)) (U^x) - g(x)g\x + ft)) f 



4V x ,fi l 
^£2Tr[/-^(x)[/ M (x)£(x + /i) 



Writing 

9 U^(x) = g\x)U^(x)g{x + ft) 
it is is clear that minimizing (|2.1|) is equivalent to maximizing 

X,fJ, 



[2.1) 



(2.2) 



(2.3) 



in the fundamental representation, which is just the lattice Landau gauge. 

Generalizing the above idea slightly, suppose we are interested in finding the best 
fit of U^(x) by a thin center vortex configuration, which has the form 



V^x) = g(x)Z fl (x)g\x + fi) 



(2.4) 



where the Z^x) = ±1 are Z2 link variables. Since the adjoint representation is blind 
to center elements, we may proceed in the following way: First find the vacuum 
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configuration g(x)g jf (x + p) which is a best fit to U^(x) in the adjoint representation, 
i.e. which minimizes 



4 = ^jE Tr A \(u,(x) - g{x)g\x + /})) (t^(x) - g{x)g\x + /})) 1 



— E 2Tr A [/ - ^(x)C^(x)^(x + jl) 

x,ix 



It is clear that the minimum of d\ is obtained at the maximum of 

R = Y.^a[ 9 U»{x)\ 

x,n 



(2.5) 



(2.6) 



which is, by definition, direct maximal center gauge. Having determined g(x) in this 
way (up to a residual Z 2 symmetry), we then want to find the Z^(x) such that the 
thin vortex configuration V^x) is a best fit to U^x). This means minimizing, in the 
fundamental representation 



d 2 = ^ E Tr \{u,(x) - g{x)Z^{x)g\x + ft)) (u^x) - g(x)Z^x)g\x + /i)) 1 



4V 



E 2Tr [ J - Z,(x)g\x)U,(x)g(x + /}) 



For fixed g(x), minimization is achieved by setting 

Z,{x) = signTrP^(x)] 



(2.7) 



(2.8) 



which is the center projection prescription. In this way we see that maximal center 
gauge, plus center projection, is equivalent to finding a "best fit" of the given lattice 
gauge field U^(x) by a thin vortex configuration V^(x). This point was first made by 
Engelhardt and Reinhardt [3] in the context of continuum Yang-Mills theory; here 
we have transcribed their argument to the lattice. 

Unfortunately, it is clear that V^(x) must be a very bad fit to U^x) at links 
belonging to thin vortices (i.e. to the P-plaquettes formed from Z^(x)). Again, this 
point was made in ref. [3] in the context of the continuum theory, on the grounds 
that the field strength of a thin vortex is divergent at the vortex core. We recall 
that a plaquette p is a P-plaquette iff Z(p) = —1 (where Z(C) denotes the product 
of Zn(x) around the contour C) and that P-plaquettes belong to P-vortices. This 
means also that 

-1 (2.9) 



^Tr[V(p)} = Z{p) 



On the other hand, at large (5 we generally have 



-Tr[U(p)] 



(2.10) 
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even at P-plaquettes. Writing link variables U^(x) in the form 



U^x) =g(x)Z ll (x)e iA ^g\x + il) 
9 U„{x) = Z^x)e lA ^ 



(2.11) 



we have 



Tr[U(p)} = Z(p)Tr[ ]J 




(2.12) 



linksdp 



or, in view of (l2~9j),(^T0l), 



\n n 




(2.13) 



But, from eq. Q2.8|) , 



Tr \ e iA ^ x) ] > 



(2.14) 



on every link. Taken together, the last two equations imply that for some links (at 
least one) on a P-plaquette, 9 U^(x) deviates strongly away from the center elements 
±1, and, as a result, 



It follows that the trace of links in the adjoint representation, which is the quantity 
being maximized in direct maximal center gauge, will be much smaller, on average, 
in the vicinity of a P-plaquette than on the rest of the lattice. 

If we compare the fits to a thermalized lattice that can be obtained from a vacuum 
configuration (no P-vortices), and from a configuration containing some arrangement 
of thin vortices on vacuum background, the latter is likely to be a better fit in the 
lattice region exterior to P-plaquettes, but a much worse fit at the P-plaquettes 
themselves. In the (3 — > oo limit, the bad fit near the P-plaquettes may overwhelm 
the good fit in the exterior region, particularly if the thermalized lattice contains 
thick 1 fm) center vortices which overlap substantially. In that case, the best fit 
is just a pure-gauge, and the global maximum of R would reveal no vortices at all. 

In our opinion, this is the explanation for the discouraging result found by BKP 
in ref. [1]. Using the simulated annealing method described in that reference, we 
have found that center projection in direct maximal center gauge has good center 
dominance properties at strong couplings, but that center dominance degrades as we 
enter the scaling regime. We have also found that Gribov copies obtained by our 
original method of over-relaxation arrive at a better fit to the thermalized lattice, 
as compared to copies generated by simulated annealing, if the computation of R 
is restricted to plaquettes in the region exterior to P-plaquettes. This is consistent 
with the idea that the bad fit in the neighborhood of P-plaquettes is the source of the 
trouble. For details and numerical results regarding this point, we refer the reader 
to ref. [4]. 




(2.15) 
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Having arrived at a diagnosis of the problem, the question is what to do about 
it. One possibility is simply to continue using the original over-relaxation technique, 
on the grounds that it provides a better fit to the lattice in the region exterior to 
P-plaquettes. Unfortunately, if one follows this idea further, gauge-fixing a large 
number of random copies on the gauge orbit and selecting the copy with the best 
fit in the exterior region, the string tension comes out too high [4]. Therefore, we 
cannot recommend the "business-as-usual" approach with much enthusiasm. 

Another possibility, suggested by Engelhardt and Reinhardt [3], is to modify R 
in eq. ( |1 . 1| ) by introducing a form factor F into the expression, i.e. 



in an effort to soften the bad fit at P-plaquettes. We have experimented with this 
form of F, but the whole approach is obviously plagued by arbitrariness. With an 
extra free parameter such as A (that can be reset at each /3), it is not surprising that 
one can adjust the projected string tension to the desired result. 

Finally, there is the Laplacian center gauge [6] developed by de Forcrand and 
Pepe [7] , which is free of Gribov ambiguities altogether. Our present article is inspired 
in large part by the de Forcrand-Pepe approach (as well as by the earlier work of 
Vink and Wiese [5]), but we do have some reservations about the version of Laplacian 
center gauge discussed in refs. [6,7]. These reservations concern the lack of scaling 
of the vortex density, the lack of precocious linearity in the vortex potential, and the 
prescription for locating vortices from the co-linearity of two Laplacian eigenvectors, 
which appears to lack the important "vortex- finding property" for thin vortices [8]. 
This last point, which is relevant to the program of looking for gauge- fixing ambi- 
guities in Laplacian gauges, will be discussed in section 5. In the meantime, we will 
proceed to the proposal which is main point of this article. 

3. Direct Laplacian Center Gauge 

To motivate our proposal, we return temporarily to the idea that, since the fit is 
worst at P-plaquettes, it might be sensible to exclude those contributions from the 
quality-of-fit function R. Specifically, consider introducing a configuration-dependent 
weighting factor 




(2.16) 



In particular, one could try the "lower bound" function 




x > A 
x < A 



(2.17) 




x G P-vortex 
otherwise 



(3.1) 
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where center projection of 9 U fl is used to decide whether or not the site x belongs to 
a P-plaquette, and c is a constant determined by the constraint 

kp^;^) = i (3.2) 

* X 

which means that c 2 is inversely proportional to the lattice volume exterior to P- 
plaquettes. One then considers 

R' = J2 P(x; 9 U,)p(x + fi- m,)TT A y(x)U,(x)g(x + p)] (3.3) 

X,fi 

which is proportional to the average trace of adjoint links Tr^U^x)] in the exterior 
volume, choosing g(x) to maximize this quantity. 

However, even the exclusion of all links joined to P-plaquettes is still not good 
enough in the continuum (3 —>■ oo limit, since there one expects a bad fit also in 
a finite volume surrounding the P-plaquettes. At low or intermediate values of /3, 
on the other hand, excluding 100% of the P-plaquette contributions goes too far, 
because the excluded links are a very substantial fraction of the total lattice volume. 
The next step, then, is to allow greater flexibility in the weighting factor. This is 
done by allowing p(x) to be a degree of freedom in its own right. Introduce a 3 x 3 
matrix-valued field 

Gij(x) = ±n ij (x)p(x) (3.4) 
where Q is an SO(3) matrix, and consider choosing p,Q to maximize the quantity 

R" = Tt[G T (x)U A „(x)G(x + p)} (3.5) 

X,fl 

where Ua/i{x) denotes the link variable in the adjoint representation, and where p(x) 
is a real positive scalar, subject to the constraint 

^EP 2 W = 1 (3-6) 

This constraint is now simply the requirement that G(x) is orthogonal "on average." 
Mapping the SO(3) matrix field fl(x) onto an SU(2) matrix field g(x) with, e.g., 
Ti[g] > 0, and again defining the lattice field in this new gauge by 9 U )M (x), one can 
determine the P- vortex locations from the center projected configuration (|2.8|) . 

Generalizing this stategy one step further, we can replace the positive scalar field 
p(x) by a real, symmetric matrix Pij(x) with positive semidefinite eigenvalues, and 
write 

M^x) = ±n ik (x)P kj (x) (3.7) 

(summation is over repeated indices). This is known as the polar decomposition of 
the matrix M. A corollary of the Singular Value Decomposition Theorem [9] is that 
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any real 3x3 matrix M can be decomposed in this way. The idea is then to find 
the 3x3 matrix-valued field M{x) which maximizes 



R M = Y / Tr[M T (x)U A ^x)M(x + 



(3.f 



X,fl 



subject to the constraint (which generalizes ( ftp] 

1 



V 



J2M t {x)M{x) 



(3.9) 



so that M(x) is also orthogonal "on average." 

It is fortunate that the problem we have just posed has both a unique solution, 
and a standard computational algorithm for arriving at that solution. It is convenient 
to view the columns of M(x) at any site set of three 3-vectors 



f b a (x) = M ab (x) , f\x) 



/$(*) 



(3.10) 



Then in this notation, maximizing Rm is equivalent to maximizing 
S = J2 f<(x)[U Alx (x)] jk f c k (x + //) + l -A cd ^ [ff{x)f}(x) - 5, 



(3.11) 



where the A is a real symmetric matrix of Lagrange multipliers, introduced to enforce 
the constraint ( |3.9j ). In view of the constraint, the configuration maximizing S also 
maximizes the expression 



S" = E [ff(x)[U A »(x)] jk f c k (x + p) - f?{x)f? 



x , 



U cd J2[f^x)ff(x)-5 cd 

x 

(3.12) 



Variation of S' wrt /f then leads to a lattice Laplacian equation 

J2^(x,y)f!(y) = AaJ7(x) 



(3.13) 



where 



T> i:j (x, y) = -^2 {pA^x^ijdy^jx + [U Afl {x - ftySy^a - 2Sij5 xy ) . (3.14) 
Because A is real and symmetric, it can be diagonalized by an orthogonal matrix O, 



i.e. 



A = T A D , Ad = diagfAi, A 2 , A 3 
JT(x) = O ab f a (x) 



(3.15) 
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and then the equation satisfied by each f' a (x) is the Laplacian eigenvalue equation 

J2^ij(x,y)f- a (y) = Kfl a (x) (no sum over a) (3.16) 
y 

Since V{j (x, y) is hermitian, the orthogonality constraint (|3.9| ) is satisfied by 
choosing the A a to be three different eigenvalues. Substituting the solutions of ( |3.16| ) 
back into (|3.12 ), its not hard to see that S' is maximized by choosing eigenvectors 



corresponding to the three lowest eigenvalues. An efficient numerical algorithm for 
obtaining the low-lying eigenvectors of a large, sparse matrix (such as our lattice 
Laplacian) is the Arnoldi method [9]. Conveniently, Fortran routines implementing 
this method, and easily adaptable to the problem at hand, are freely available [10]. 

Since 5" is invariant under the global transformation — > f- a (x) = O a bf^(x), 
we may as well use the solution for the maximum for which M a b(x) = fa(x). At this 
point, we have to map the matrix field M(x) onto an SO(3)-valued field 5u(x), and 
this in turn to a corresponding SU(2) gauge transformation g(x). We consider two 
possibilities: 

1. Naive Map: Choose Qa{%) to be the SO(3)-valued field which is closest to 
M(x), in the sense that 

Tr[^(x)M T (x)]| (3.17) 

is maximized at each site x. 

2. Laplacian Map: Choose [^(^)]y — fi( x ) to be the SO(3)-valued field closest 
to M(x) subject to the constraint that #a(x), like M(x), satisfies a Laplacian 
equation of the form 

= ^ac(x)mx) (3.i8) 

y 

The "naive" choice is equivalent to the Laplacian Landau gauge of Vink and 
Wiese [5], generalized to the adjoint representation. It is obtained by first making a 
singular value decomposition 

M(x) = U(x)M D (x)V T (x) (3.19) 

with U, V orthogonal matrices and Mp a diagonal positive semi-definite matrix. 1 
Then make the polar decomposition at each site 

M{x) = [U(x)V T (x)}[V(x)M D (x)V T (x)} = ±g A (x)P(x) (3.20) 

with gA.{x) the SO (3) matrix- valued field 

g A {x) = det[U(x)V T (x)]U(x)V T (x) (3.21) 

^^Again, there are standard numerical packages which implement singular value decomposition. 
We have used the routine svdcmp from Numerical Recipes [11]. 
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and P is symmetric and positive semi-definite. Matrix g A (x) is guaranteed to be the 
S0(3) matrix closest to M(x), in the sense of maximizing the expression ( |3.17|) [9]. 



Finally, we map qa^x) at each site onto one of its two possible representatives in the 
SU(2) group 



I + a a a b [g A (x)} 



g(x) = ± \' — f Av r y A°f (3.22) 
yy } 2[l + Tr(g A (x))] 1 ' 



The choice of sign in ( |3.22| ) is irrelevant, all choices being related by the residual 
Z 2 invariance. One can then transform the lattice configuration by g(x) to obtain 
9 Ufj,(x), and carry out center projection. 

The drawback of the "naive" choice is that, while M(x) is a covariantly smooth 
matrix, whose columns are the low-lying eigenvectors of the lattice Laplacian, this 
covariant smoothness property is less pronounced in the gA{x) obtained from singular 
value decomposition. This results in high-frequency "noise" (small scale fluctuations) 
in the P-vortex surfaces of the projected configuration. We have found the Laplacian 
mapping to be preferable. 

To carry out the Laplacian mapping, we begin with the "naive" mapping, to find 
the SO(3) matrix field Q(x) closest to the real- valued matrix field M(x) obtained from 
solving the Laplacian eigenvalue equation fl3.l6| ). Then we relax fl(x) to the closest 
solution 5u(x) satisfying the generalized Laplacian equation ( 3.18| ). The relaxation 



is carried out in the following way: We note that the equation (|3.18|) is simply the 
stationarity condition of the action 

S" = E [f-(x)[U All (x)} jk r k (x + p)- tf(x)ff{x)] +W A cd (x) {./J(.r)/;'(.r) - 5 cd } 

(3.23) 

with [gA(x)]ij = fi( x ) an SO(3) field. With the local Lagrange multiplier field 
enforcing the SO(3) constraint, a stationary solution of 5"' is easily seen to be a 
Gribov copy of direct maximal center gauge; local maxima of S" are also local maxima 
of R in eq. Q. 

If S[g] = —R is regarded as the action of a spin system, with g A {x) the SO(3)- 
valued field variables, then each Gribov copy can be regarded as a metastable state 
of system, each with its own "basin of attraction." The basin of attraction of a 
metastable state is the volume of all configurations which, when the system is sud- 
denly cooled (or "quenched"), will fall into that metastable state. Over-relaxation is 
essentially a sudden cooling technique [4] , and is therefore perfectly suited to sliding 
the system "down the hill" of action S, from the configuration Q UAfi{x) to the nearest 
(or, at least, a nearby) local minimum of S (see Fig. I). 2 The gauge transformation 
g(x) obtained at this minimum is the Laplacian map of M(x), and the corresponding 
9 U tl {x) is the lattice configuration in direct Laplacian center gauge. 



2 There is no guarantee that the local minimum obtained in this way is truly the nearest of the 
local minima to M(x), as this depends on the topography of S[g] in the neighborhood of f2. 



10 



PD ... Polar Decomposition 
OR ... Over-Relaxation 



Figure 1: Laplacian mapping of M{x) to a nearby S0(3) matrix-valued field gA(%)- 

It is interesting that a drawback of the over-relaxation technique, in fixing to 
maximal center gauge, becomes a virtue in direct Laplacian center gauge. The draw- 
back is that over-relaxation is not good at finding the global minimum of S, but 
only takes a given initial state to a nearby local minimum (Gribov copy). In fixing 
to direct Laplacian gauge, on the other hand, we have taken advantage of the fact 
that over-relaxation goes to a nearby minimum in order to carry out the Laplacian 
mapping M(x) —> g(x). 

It may be useful at this stage to summarize the steps of direct Laplacian center 
gauge fixing: 

1. From a thermalized SU(2) lattice, construct the lattice configuration in the 
adjoint representation 

[U A ,(x)] l3 = hrrfaU^xtoUlix)] (3.24) 



2. Solve the eigenvalue problem (|3.16 ) for the eigenvectors corresponding to the 



three smallest eigenvalues of the lattice Laplacian operator in the adjoint rep- 
resentation. We have used the ARPACK routines [10] for this purpose. From 
the eigenvectors f a (x), a = 1,2,3, construct the matrix field M^-(x) = fl{x). 

3. Perform, at each site, the singular value decomposition M(x) = U(x)Md(x)V t (x) 
and extract the SO (3) matrix- valued field 

Q(x) = det[U(x)V T (x)}U(x)V T (x) (3.25) 

4. Map Q to an SU(2) matrix- valued field 

, , I + G a a b \tt(x)\ ab 
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and perform the gauge transformation 

"U^x) = J{x)U^(x)uj(x + p) (3.27) 

The result is to transform the configuration into Laplacian Landau gauge in 
the adjoint representation. 

5. With "U^x) as the starting point, relax the lattice field to the nearest (or, at 
least, a nearby) Gribov copy of direct maximal center gauge, using the over- 
relaxation method described in ref. [2]. The result is the configuration fixed to 
direct Laplacian center gauge. 

To make a long story short, direct Laplacian center gauge in practice is nothing 
but gauge fixing to adjoint Laplacian Landau gauge, followed by over-relaxation to 
the nearest Gribov copy of direct maximal center gauge. Both the direct Laplacian 
and the direct maximal center gauges aim at generating a configuration satisfying 
locally the adjoint lattice Landau gauge condition. The difference is that the opti- 
mal configuration, in maximal center gauge, is generated by a gauge transformation 
maximizing the R functional, whereas in direct Laplacian center gauge the optimal 
configuration is generated by a gauge transformation lying closest to a matrix maxi- 
mizing the Rm functional. We hope to have explained, in the previous section, why 
the global maximum of R is not necessarily the best choice for vortex finding, and to 
have motivated the alternative choice of Gribov copy corresponding to direct Lapla- 
cian center gauge. In the next section, we present results that are obtained from 
center projection in this gauge. 

4. Numerical Results 

To test the reasoning of the last section, we have recalculated the vortex observables 
introduced in our previous work (cf. ref. [2,12]), with P-vortices located via center 
projection after fixing the lattice to the new direct Laplacian center gauge. 

4.1 Center Dominance and Precocious Linearity 

In view of the ~ 30% breakdown of center dominance in maximal center gauge found 
by BKP [1], the quantities of most immediate interest are the center-projected Creutz 
ratios, extracted from I x J Wilson loops on the center-projected lattice. These are 
displayed, for (3 in the range [0.4,2.5], in Fig. ^|. 

At strong couplings, it is clear from Fig. |2| that Xcp(l ; 1) matches up with the 
analytic prediction, in the full theory, of a = — ln(/5/4). At weaker couplings, we 
display our data for (3 = 2.4 and (3 = 2.5 in Figs. |3| and |j. Once again we see the 
feature of precocious linearity; i.e. the fact that Creutz ratios Xcp(R,R) at a given 
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Figure 2: Creutz ratios from center-projected lattice configurations, in the direct Lapla- 
cian center gauge. 

(3 vary only slightly with R, which means that the center-projected potential is ap- 
proximately linear starting at R = 2 lattice spacings. The significance of precocious 
linearity is that it implies that the center-projected degrees of freedom have isolated 
the long-range physics, and are not mixed up with ultraviolet fluctuations. 

It is also apparent that the center-projected Creutz ratios are quite close to the 
asymptotic string tension of the unprojected theory, reported in refs. [13]. This is 
evidence, at least at these two couplings, of the center dominance property. Our 
combined data for the range of couplings j3 = 2.2 — 2.5 is displayed on a logarithmic 
plot in Fig. [5|. In general Xc P {R, R) deviates from the full asymptotic string tension 
by less than 10%. 

As another way of displaying both center dominance and precocious linearity, we 
adopt the usual procedure of assigning a lattice spacing a(/3) based on the asymptotic 
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(3=2.4 



0.12 




Figure 3: Center projected Creutz ratios at f3 = 2.4 on various lattice sizes. The horizontal 
band indicates the full asymptotic string tension and errorbar. 
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Figure 4: Same as Fig. |j| at (3 = 2.5. 



string tension in lattice units o" La< (/5), and the string tension in physical units a p h ys 
(440 MeV) 2 , i.e. 

<?Lat(f3) 



&phys 



(4.1) 
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Figure 5: Combined data, at (3 = 2.2 — 2.5, for center-projected Creutz ratios obtained 
after direct Laplacian center gauge fixing. Horizontal bands indicate the asymptotic string 
tensions on the unprojected lattice, with the corresponding errorbars. 

We then display, in Fig. || the ratio 

x P hys(R, R) _ Xc P (R, R) ^ 2 ^ 

&phys &Lat 

as a function of the distance in physical units 

R phys = Ra{(3) (4.3) 

for all Xc P (R, R) data points taken in the range of couplings (3 = 2.3 — 2.5. Again 
we see that the center-projected Creutz ratios and asymptotic string tension are in 
good agreement (deviation < 10%), and there is very little variation in the Creutz 
ratios with distance. 
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Figure 6: The ratio of projected Creutz ratios to the full asymptotic string tension, as 
a function of loop extension in fermis. The data is taken from Xcp(R,R) at a variety of 
couplings and lattice sizes. 



Precocious linearity can be understood in the following way: Asymptotically, 
outside the vortex core, a center vortex is a Zn dislocation, which affects the value 
of large Wilson loops, topologically linked to the vortex, by a center element factor. 
Assume that P-vortices correctly identify the middle of thick center vortices in lattice 
configurations. Then center projection in effect collapses the thick (ps one fermi) core 
of a vortex to a width of one lattice spacing. This means that the asymptotic effect 
of thick vortices is obtained, on the projected lattice, at any distance greater than 
one lattice spacing. If P-plaquettes in a plane are uncorrelated, this leads to a 
linear potential at short distances on the projected lattice. On the other hand, if 
precocious linearity is not found on the projected lattice, it means that either the P- 
vortex surface is very rough, fluctuating on all distance scales, or else that some large 
fraction of P-plaquettes on the projected lattice belong to P-vortices which are small 
in extent, and do not percolate. Either case results in some short-range correlations 
among P-plaquettes in a plane, corresponding to high-frequency phenomena not 
directly associated with the long-range physics. 



4.2 Scaling of the Vortex Density 

Denote the total number of plaquettes on the lattice by N T) the number of P- 
plaquettes by N vor , and the density N vot /Nt of P-plaquettes on the lattice by p. 
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We have 



-t* vor ^vor^ 1 2 

~Nt ~ iV T a 4 a 
Total Vortex Area 
6 x Total Volume 




where p is the center vortex density in physical units. Then, according to asymptotic 
freedom 
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6A 2 




67T 
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(4.4) 



The average vortex density p can be extracted from the expectation value of center 
projected plaquettes 

W cp (l, 1) = (1 - p) x (+1) +px (-1) = l-2p 

p= l -{l-W cp {l ) l)) (4.5) 

Figure [7] is a logarithmic plot of P-vortex density p vs. (3. Errorbars are less 
than the size of data points. The solid line is the asymptotic freedom prediction 



with J p/QA 2 = 50. We emphasize that the slope of this line represents the 
proper asymptotic scaling for surface densitites. The slope that would be associated 
with the scaling of pointlike objects such as instantons, or linelike objects such as 
monopoles, would be quite different. Note that the apparent scaling of Xcp(l, 1) = 
— In W cp (l, 1) ~ 2p in Fig. || is just a consequence of the scaling of p. 

4.3 Vortex-Limited Wilson Loops 

A "vortex- limited" Wilson loop W n {C) is defined as the expectation value of an 
unprojected Wilson loop around some contour C, evaluated in the sub-ensemble of 
configurations in which, on the corresponding center-projected lattice, precisely n 
P-vortices pierce the minimal area of the loop. If P-vortices on the projected lattice 
roughly locate the middle of thick center vortices on the unprojected lattice, then in 
the limit of large loop areas, for the SU(2) gauge group, we expect [2, 12] 

The numerical evidence definitely shows a trend in this direction, as can be seen from 
our data for W\/Wq and W^/Wq vs. minimal loop area at (3 = 2.3, shown in Figs. || 
and |9|. 
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p 

Figure 7: Evidence of asymptotic scaling of the P-vortex surface density. The solid line 
is the asymptotic freedom prediction (p~4|), with y/ p/6A 2 = 50. 

If P-vortices in the projected configuration locate thick center vortices in the 
unprojected configurations, and if center vortices are responsible for confinement, 
then Wq(C) should not have an area-law falloff. In Fig. [10] we compare Creutz ratios 
Xo[I, J] extracted from rectangular Wo(I, J) loops, with the standard Creutz ratios 
x(I, J) from loops evaluated in the full ensemble. As expected, the Creutz ratios of 
the zero-vortex loops tend to zero with loop area. 

4.4 Center Vortex Removal 

It was suggested by de Forcrand and D'Elia [14] that one could remove center vortices 
from a given lattice configuration by simply multiplying that configuration by the 
corresponding center-projected configuration derived in maximal center gauge, i.e. 

Ufa) = Zfa)Ufa) (4.7) 
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16 4 Lattice at (3 = 2.3 




-0.5 



10 15 
Loop Area 



Figure 8: Ratio of the one- vortex to zero- vortex Wilson loops W\{C) /Wq(C) vs. loop 
area, at /3 = 2.3 on a 16 4 lattice. 
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Figure 9: Ratio of the two- vortex to zero- vortex Wilson loops W2(C)/Wq(C) vs. loop 
area, at j3 = 2.3 on a 16 4 lattice. 



where Z^(x) is given by ( |2.8|) . Since the adjoint representation is blind to center 
elements, it is easy to see that if g(x) is a transformation such that 9 U fJi (x) is in 
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Figure 10: Creutz ratios xo(RjR) extracted from zero-vortex Wilson loops Wo(I, J), as 
compared to the usual Creutz ratios x(R, R) on a 16 4 lattice at = 2.3. 

maximal center gauge, then 9 U' (x) is also in maximal center gauge. However, there 
are no P-vortices obtained from center projection of 9 U' since 

%{x) = affOt[%(x)] = Z\{x) = 1 (4.8) 

One can therefore say that center vortices, as identified in maximal center gauge, have 
been removed from the lattice configuration. More precisely, what the modification 
( f4.7p does is to place a thin vortex (one plaquette thickness) in the middle of each 
thick center vortex core, whose locations are identified by center projection. At large 
scales, the effects of the thin and thick vortices on Wilson loops will cancel out. Thus 
there should be no area law due to vortices in the modified configuration U' and 
the asymptotic string tension should vanish. 

The vanishing of string tension in the modified configurations was, in fact, ob- 
served in ref. [14], using maximal center gauge fixing by the over-relaxation technique 
of ref. [2]. Figure [11] is a repeat of the de Forcrand-D'Elia calculation at (3 = 2.3, 
using direct Laplacian center gauge rather than direct maximal center gauge, and 
we find essentially the same result. 

5. Remarks on Gauge-Fixing Ambiguities 

The original suggestion of 't Hooft, in his 1981 article on monopole confinement [15], 
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Figure 11: Creutz ratios on on the modified lattice, with vortices removed, at (3 = 2.3. 
For comparison, we also display the unprojected Creutz ratios (open circles), the center 
projected Creutz ratios (solid circles), and the asymptotic string tension (horizontal band). 

was to introduce a composite operator in the pure gauge theory transforming like 
a field in the adjoint representation of the gauge group; for SU(2) gauge theory 
this operator defines a unitary gauge leaving a residual U(l) symmetry. Monopole 
worldlines would then be associated with lines along which the gauge transformation 
is ambiguous. This was the idea which motivated the numerical study of maxi- 
mal abelian [16] and Laplacian abelian [17] gauges. The idea can be generalized by 
introducing not one but two composite operators transforming in the adjoint repre- 
sentation. Let us denote these operators, in SU(2) gauge theory, by 4> a (x) and r} a (x), 
with color index a = 1, 2, 3. One then performs a gauge transformation which takes, 
e.g., <p a into the positive color 3 direction at every point, and r\ a to lie in the color 
1-3 plane, with 03, 771 > 0. This gauge leaves a remnant Z<i symmetry. Gauge-fixing 
ambiguities occur on surfaces where (j>(x) and r](x) are co-linear, and these surfaces 
are then to be identified with center vortices. 

The suggestion of refs. [6,7] is to choose, for 4> a (x) and r) a (x), the two lowest 
lying eigenstates of the covariant lattice Laplacian operator ( |3.14| ) in the adjoint 
representation; i.e. 

= , Va(x) = f a (x) (5.1) 

with Ai < A2 being the two lowest eigenvalues of the Laplacian operator. This is the 
original version of Laplacian center gauge, which we will refer to as LCG1. Because 
the procedure involves first fixing to Laplacian abelian gauge, followed by a further 
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gauge fixing which reduces the residual gauge symmetry from U(l) to Z 2 , LCG1 is 
reminiscent of the indirect maximal center gauge of ref . [12] . As we have pointed out 
in refs. [19], abelian monopole worldlines in the indirect maximal center gauge lie 
on center vortex surfaces, and a vortex at fixed time can be viewed as a monopole- 
antimonopole chain. The same relationship between abelian monopoles and center 
vortices holds true in LCG1. 

Gauge-fixing ambiguities in LCG1 occur when and 7] are co-linear in color 
space, and it is suggested that these ambiguities can be used to locate vortex sur- 
faces. This approach to vortex finding is certainly quite different from the reasoning 
outlined in section 3, which is is motivated by a "best fit" procedure. In special 
cases, notably for a classical vortex solution on an asymmetric lattice with twisted 
boundary conditions, the co-linearity approach seems to work well [18]. We take 
note, however, of a simple counter-example: Suppose we insert two vortex sheets 
"by hand" into a configuration U^{x) by the transformation 



All other components are unchanged, i.e. U' k (x) = Uk{x) for k > 0. This transforma- 
tion inserts two thin vortex sheets into the lattice, parallel to the x<i — x 3 plane, at 
xo = 0. We then ask whether these two vortex sheets will be located by the gauge- 
fixing ambiguity approach. The answer is clearly no, since the composite Higgs fields 



urations. This example illustrates the fact that the gauge-ambiguity approach lacks 
the "vortex-finding property" discussed in ref. [8] , which is the ability of a procedure 
to locate thin vortices inserted at known locations into the lattice. We have argued 
in ref. [8] that this property should be a necessary (although not sufficient) condition 
for locating vortices on thermalized lattices. 

In practice, on thermalized lattices, vortices are located in LCG1 via center 
projection, rather than by eigenvector co-linearity [7]. Using center projection to 
locate vortices, LCG1 recovers the vortex-finding property, for reasons explained in 
ref. [8]. 3 It also exhibits center dominance of the projected asymptotic string tension. 
On the other hand, the projected Creutz ratios in LCG1 do not display precocious 
linearity [7], nor does the vortex surface density scale according to the asymptotic 
freedom formula [20]. A possible remedy is to first fix the lattice to LCG1, and 
from there to fix to (direct or indirect) maximal center gauge by an over-relaxation 
procedure. This latter procedure has, in fact, been tried by Langfeld et al. in ref. [20], 
with good results for the vortex density. Other vortex observables, discussed in the 
previous section, have not yet been studied systematically in this approach, which 

3 As discussed in that reference, the vortex-finding property is obtained from center projection 
in any gauge which completely fixes link variables in the adjoint representation, leaving a residual 
Zn invariance. 




xq = and < x\ < L 
otherwise 



(5.2) 



(f)(x) and i](x) are identical in the original U^(x), and in the modified U' (x) config- 
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seems to have a great deal in common with the direct Laplacian center gauge we 
have advocated here. 



6. Conclusions 

We have tested a procedure for locating center vortices on thermalized lattices, based 
on the idea of finding the best fit to the thermalized lattice by thin vortex configu- 
rations. Our new procedure, which is essentially just a variation of direct maximal 
center gauge, is designed to soften the inevitable bad fit at vortex cores, due to the 
singular field strength of thin vortices. The numerical results we have found are 
promising: Deviations from center dominance are generally less than 10%, P-vortex 
density scales correctly, and there are the usual strong correlations between P-vortex 
locations and gauge-invariant observables. Our new method is motivated by an im- 
proved understanding of how maximal center gauge works, and it addresses some 
objections to center gauge fixing that have been raised in the recent literature. We 
hope it will provide a more solid foundation for further numerical studies. 
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